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Abstract — This paper considers cooperative spectrum sensing 
in Cognitive Radios. In our previous work we have developed 
DualSPRT, a distributed algorithm for cooperative spectrum 
sensing using Sequential Probability Ratio Test (SPRT) at the 
Cognitive Radios as well as at the fusion center. This algorithm 
works well, but is not optimal. In this paper we propose an 
improved algorithm- SPRT-CSPRT, which is motivated from Cu- 
mulative Sum Procedures (CUSUM). We analyse it theoretically. 
We also modify this algorithm to handle uncertainties in SNR's 
and fading. 

Keywords- Cognitive Radio, Spectrum Sensing, Cooperative 
Distributed Algorithm, SPRT. 

I. Introduction 

Presently there is a scarcity of wireless spectrum worldwide 
due to an increase in wireless services. Cognitive Radios 
are proposed as a solution to this problem. They access the 
spectrum licensed to existing communication services (pri- 
mary users) opportunistically and dynamically without causing 
much interference to the primary users. This is made possible 
via spectrum sensing by the Cognitive Radios (secondary 
users), to gain knowledge about the spectrum usage by the 
primary devices. However due to the strict spectrum sensing 
requirements [21] and the various wireless channel impair- 
ments spectrum sensing has become the main challenge faced 
by the Cognitive Radios. 

Cooperative spectrum sensing ([24], ED) in which different 
Cognitive Radios communicate each other exploits spatial 
diversity among them effectively. This can largely solve the 
problems caused by shadowing, multipath fading and hidden 
node problem in spectrum sensing. Moreover it improves the 
probability of miss detection and the probability of false alarm. 
Cooperative spectrum sensing ( 071 . ll27ll ) is called centralized, 
when a central unit gathers sensing data from the Cognitive 
Radios and identifies the spectrum usage. It is distributed if 
each local user uses the observations to make a local decision 
and sends this to the fusion center to make the final decision. 
Secondary users can either transmit a soft decision (summary 
statistic) or a hard decision [ 17|. Soft decisions provide better 
performance but at the cost of higher bandwidth consumption 
by the control channels between the Cognitive Radio and the 
fusion center. However as the number of cooperative users 
increases, hard decisions can perform as well 0. 
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An extensive survey of spectrum sensing methods is pro- 
vided in [27|. Recent spectrum sensing schemes are based 
on higher order statistics fPfl . wavelet transforms [22] and 
compressed sensing [23]. One can use a fixed sample size (one 
shot) detectors or sequential detectors (J6), ED, ED, ll27l ). 
In fixed sample size detectors, the matched filter is optimal 
when there is complete knowledge of the primary signal. When 
the only a prior knowledge is about the noise power, then 
an energy detector is optimal in Neyman-Pearson criterion 
0. However sequential detectors perform better. A recent 
survey is |fl2l . The sequential detectors can detect change 
or test a hypothesis. Sequential hypothesis testing finds out 
whether the primary is ON or OFF, while the sequential change 
point detection detects the point when the primary turns ON 
(or OFF). Sequential change detection is well studied (see 
10, ED, 021, 03 and the references therein). However the 
optimal solution in the distributed setup is still not available. 
Sequential hypothesis testing (0, 0, ED, ED) is useful 
when the status of the primary channel is known to change 
very slowly, e.g., detecting occupancy of a TV channel. Usage 
of idle TV bands by the Cognitive network is being targeted as 
the first application for cognitive radio. In this setup Walds' 
Sequential Probability Ratio Test (SPRT) [25 1 provides the 
optimal performance for a single node ( fl6l . [26]). But the 
cooperative setup is not well explored. 

We consider cooperative spectrum sensing using sequential 
hypothesis testing. SPRT is used at both the secondary nodes 
and the fusion center. This has been motivated by our pre- 
vious algorithm DualCUSUM for change detection [6|. This 
algorithm is called DualSPRT and has been studied in (Q, 
and ED). As against and ED, in it has been 
analysed theoretically also and has been extended to cover 
channel and noise power uncertainties. Cooperative spectrum 
sensing via sequential detection is also considered in [28 1. 
But fusion center noise is not modelled in 1281 . Similarly 
03l provides the optimal decentralized sequential hypothesis 
testing algorithms without considering fusion center noise. 
Neither does it consider SNR uncertainty and fading. 

Although DualSPRT works well, it is not optimal. In this 
paper we improve over DualSPRT. Furthermore we introduce 
a new way of quantising the local nodes SPRT decisions. We 
call this algorithm SPRT-CSPRT. We extend this algorithm to 
cover SNR uncertainties and fading channels. We also provide 



its theoretical analysis. 

This paper is organised as follows. Section II describes 
the model. Section III starts with the DualSPRT algorithm. 
Then we provide SPRT-CSPRT and DualCSPRT algorithms 
developed in this paper. We compare their performance. Next 
we consider the receiver SNR uncertainty and slow fading 
channels. Section IV provides the theoretical analysis and 
compares to simulations. Section V concludes the paper. 

II. Model 

Consider a Cognitive Radio system with one primary trans- 
mitter and L secondary users. The L local nodes sense the 
channel to detect the spectral holes. The decisions made by the 
secondary users are transmitted to a fusion node via a Multiple 
Access Channel (MAC) for it to make the final decision. 

Let Xk,i be the observation made at secondary user I at 
time k. We assume that {Xk,i, k > 1} are independent 
and identically distributed (i.i.d.) and that the observations 
are independent across Cognitive Radios. Using the detection 
algorithm based on {X n j, n < k} the secondary user I 
transmits Yj. ; to the fusion node. We also assume that the 
secondary nodes are synchronised so that the fusion node 
receives Yfe = Yli=i ^k,J + where {Z^} is i.i.d. receiver 
noise, it will be assumed to be zero mean Gaussian with 
variance a 2 . The fusion center observes {Yfe} and decides 
upon the hypothesis. 

The observations {Xk.i] depend on whether the primary is 
transmitting (Hypothesis Hi) or not (Hypothesis Ho): 

x _ f Z k ,i, k = 1,2,..., under H , 

k ' 1 \ hiS k + Z k ,j, fe = l,2, under H l7 

where hi is the channel gain of the I th user, Sk is the primary 
signal and Zk.i is the noise at the I th user at time k. We assume 
{Zk.i-, k > 1} are i.i.d. . Let the fusion center makes a decision 
at time N . We assume that N is much less than the coherence 
time of the channel so that the slow fading assumption is valid. 
This means that hi is random but remains constant during the 
spectrum sensing duration. 

The general problem is to develop a distributed algorithm 
in the above setup which solves the problem: 

rmnE DD =E[N\Hi] , (2) 

subject to Pfa < ct 

where Hi is the true hypothesis, i = 0,1 and Pfa is the 
probability of false alarm, i.e., probability of making a wrong 
decision. We will separately consider E[N\H{\ and E[N\Hq]. 
It is well known that for a single node case (L — 1) Wald's 
SPRT performs optimally in terms of reducing E[N\Hi] and 
E[N\Hq] for a given Pfa - Motivated by the good performance 
of DualCUSUM in (Q], 0) and the optimality of SPRT for 
a single node, we proposed DualSPRT in [7| and studied its 
performance. Now we modify DualSPRT to SPRT-CSPRT and 
DualCSPRT and we present the theoretical analysis of this 
algorithms. 



III. Sequential Spectrum Sensing Algorithms 

We first present DualSPRT which was introduced in our 
previous work J7J- 

A. DualSPRT Algorithm 

1) Secondary node, I, runs SPRT algorithm, 

W ,i = 

W kJ = W k -!,i + log [f h i (X kJ ) Ih.i (X k .i )] , k > 1 (3) 

where f±j is the density of Xk.i under Hi and/o,; is the 
density of X k> i under Hq. 

2) Secondary node I transmits a constant b\ at time k if 
Wk,i > 71 or transmits bo when Wk.i < 7o, i-e., 

Yk : i = M{iy M > 7 i} + bo^{w k ,i< 10 } 

where 70 < < 71 and 1a denotes the indicator 
function of set A. Parameters &i,&Oj7i>7o aie chosen 
appropriately. 

3) Physical layer fusion is used at the Fusion Centre, i.e., 
Yfe = Y^d=i Yk,l + %k, where Z k is the i.i.d. noise at the 
fusion node. 

4) Finally, Fusion center runs SPRT: 

Ffe = Ffe_ 1 +log[ 5l (Yfe)/ ff0 (Yfc)], F o = 0, (4) 

where go is the density of Zk + fio and gi is the density 
of Zk + Hi, Ho and /zi being design parameters. 

5) The fusion center decides about the hypothesis at time 
N where 

N = inf{fc : F k > fa or F k < /3 } 

and Po < < fa. The decision at time N is Hi if 
Fn > Pi', otherwise Ho. 

B. SPRT-CSPRT Algorithm 

In DualSPRT given above, observations to the fusion center 
are not always identically distributed. Till the first transmission 
from secondary nodes, these observations are i.i.d. ~ A/"(0, a 2 ) 
where J\f(a, b) is the Gaussian pdf with mean a and variance 
b. But after the transmission from the first local node and 
till the transmission from the second node, they are i.i.d. 
Gaussian with another mean and same variance a 2 . Thus the 
observations at the fusion center are no longer i.i.d. . Since 
the optimality of SPRT is known for i.i.d. observations ( l26ll . 
lfT6l ). DualSPRT is not optimal. 

The following heuristic arguments provide the motivation 
of the proposed modifications to DualSPRT. A sample path of 
the fusion center SPRT under the hypothesis Hi is given in 
Figure Q] If the SPRT sum defined in (0]i goes below zero it 
delays in crossing the positive threshold Pi . Hence if we keep 
SPRT sum at zero whenever it goes below zero, it reduces 
E DD . This happens in CUSUM (031, CD). Similarly one 
can use a CUSUM type algorithm under Ho. Thus we obtain 
the following algorithm, 

Steps (l)-(3) are same as in DualSPRT. The steps (4) and 
(5) are replaced by 



4) Fusion center runs two algorithms: 



Fk = (FLi + log [.9i (n) /go (Yk )] + Di)+ (5) 



*k = (Fk-i + log [5i (Y k ) /go (Y k )] + D )- , (6) 

F X = 0,F ° = 0, where (x)+ = max(0,z) and 
(x)~ = min(0,a;). D\ and Dq are appropriately chosen 
constants to introduce bias to the drift. 
5) The fusion center decides about the hypothesis at time 
N where 

N = mf{fe : Fl > fa or F% < /3 } 

and /3q < < /3i. The decision at time N is Hi if 
-Fjv > /^i> otherwise Ho. 

Under H\, © has a positive drift and hence it approaches the 
threshold j3\ quickly, but under Hq, (0 will most probably be 
hovering around zero. Similarly under Hq, ^ moves towards 
Po, but under Hi will be mostly around zero. This means that 
Pfa for this algorithm is expected to be less compared to 
DualSPRT. 
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Fig. 1. Sample Path of F k under SPRT Sum and CSPRT Sum for 71 

j3\ = 20, /ii = 1 and fio = — 1 



We consider one more improvement. When a local Cogni- 
tive Radio SPRT sum crosses its threshold, it transmits 61/60. 
This node transmits till the fusion center SPRT crosses the 
threshold. If it is not a false alarm, then its SPRT sum keeps 
on increasing (decreasing). But if it is a false alarm, then the 
sum will eventually move towards the other threshold. Hence 
instead of transmitting 61/ 60 the Cognitive Radio can transmit 
a higher / lower value in an intelligent fashion. This should 
improve the performance. Thus we modify the step (3) in 
DualSPRT as follows. Secondary node I transmits a constant 
from {6j, 62, 63, 64} at time k if W k ,i > 7i or transmits from 
{6?, b% 6°, 64*} when W k .i < 7o , as follows : 



Y k , = < 



b{ if W k j 
b\ if W k ,i 
b\ if W k ,i 
61 if W Kl 
b°i if W Kl 
b\ if W kJ 
b° 3 if W kil 
bl if W Kl 



G [71,71 +2Ai), 

G [71 +2Ai, 7 i +4Ai), 

G [71 +4Ai,7i +6A1), 

G [71 +6Ai,oo), 

G [70,70 -2A ), 

G [ 7o -2A ,7o-4Ao), 

G [70 -4A ,7o -6A ), 

G [70 - 6A , -00). 



(7) 



where Ai and Ao are the parameters to be tuned at the 
Cognitive Radio. The expected drift under Hi (Hq) is a good 
choice for Ai (Ao). 

We call the algorithm with the above two modifications as 
SPRT-CSPRT (with 'C as an indication about the motivation 
from CUSUM). 

If we use CSPRT at both the secondary nodes and the fusion 
center with the proposed quantisation methodology (we call it 
DualCSPRT) it works better as we will show via simulations 
in Section III C. In the Section IV we will theoretically analyse 
SPRT-CSPRT. As the performance of DualCSPRT (Table H is 
nearer to that of SPRT-CSPRT, we analyse only SPRT-CSPRT. 

C. Performance Comparison 

Throughout the paper we use 71 = —70 = 7, /3i = — /?o 
and fii = —fio = fi for the simplicity of the simulation and 
analysis. 

We apply DualSPRT, SPRT-CSPRT and DualCSPRT on the 
following example and compare their Ejjd for various values 
of Pfa- We assume that the pre-change distribution /q and 
the post change distribution /j are Gaussian with different 
means. This type of modelling is relevant when noise and 
interference are log-normally distributed [24|. This is also a 
useful model when X k j is the sum of energy of a large number 
of observations at the secondary node at low SNR. 

Parameters used for simulation are as follows: There are 5 
nodes (L = 5), / ,; ~ W(0, 1), for 1 < I < L. Primary to 
secondary channel gains are 0, -1.5, -2.5, -4 and -6 dB re- 
spectively (the corresponding post change means of Gaussian 
distribution with variance 1 are 1, 0.84, 0.75, 0.63 and 0.5). 
We assume Z k - J\f(0, 5) and drift of DualSPRT and SPRT- 
CSPRT at the fusion center is taken as 2fiY k , with /1 being 
1. We also take D Q = D t = 0, {6}, 6^,6^,61} = {1,2,3,4}, 
{6?, bl, bl, bl) = {-1, -2, -3, -4} and 6i=-6 =l(for Dual- 
SPRT). Parameters 7 and (3 are chosen from a range of values 
to achieve a particular Pfa- Table U provides the Ebd and 
Pfa via simulations. We see a significant improvement in 
Edd compared to DualSPRT. The difference increases as Pfa 
decreases. The performance under Hq is similar. 





Pfa = 0.1 


P FA = 0.001 


Pfa = 5e - 5 


DualSPRT 


19.74 


31.37 


34.177 


SPRT-CSPRT 


15.52 


22.59 


23.673 


DualCSPRT 


14.96 


21.52 


21.88 



TABLE I 

Comparison among DualSPRT, SPRT-CSPRT and DualCSPRT for 

DIFFERENT SNR'S BETWEEN THE PRIMARY AND THE SECONDARY USERS, 
UNDER Hi 



D. Unknown Received SNR and Fading 

In this section, now we consider the following setup. We use 
energy detector at the Cognitive Radios, i.e, the observations 
Xf.,1 are a summation of energy of past A^ observations 
received by the I th Cognitive Radio node. Then if N\ is 
reasonably large, Xk,i are approximately Gaussian. If the 
received SNR at the Cognitive Radio is not known then the 
hypothesis testing problem can be approximated as a change 
in mean of Gaussian distributions problem, where the mean 9\ 
under Hi is not known. For this case in [7| we used composite 
sequential hypothesis testing proposed in iTPTl at the secondary 
nodes and used SPRT at the fusion node. This was called GLR- 
SPRT 0. Here, to take the advantage of CSPRT at the fusion 
node and the new quantisation technique we modify GLR- 
SPRT to GLR-CSPRT with appropriate local quantisation. 
Thus the secondary node's hypothesis testing problem, SPRT, 
stopping criteria and decision are modified as follows, 



> 



W n i = max 



(8) 

(9) 
(10) 



Hq : = do ; Hi 
where 9q = and 9i is appropriately chosen, 

1^ l0g f (Y V 2-, l0g f (Y \ 

N = inf {n : W n ,i > g{cn)} , 

where g() is a time varying threshold and c is the cost 
assigned for each observation. Its approximate expression 
is given in ifTTl . Also for Gaussian fo and fx, 9 G [01,02] 
and S n as the summation of observations Xk,i upto time 
n, 9 n = max{ai,inin[S n /n,a2}}- At time N decide upon 
H or Hi according as 9n < 9* or §n > 9* , where 9* is 
obtained by solving 1(9*, 9 ) = 1(6*, 61), and 1(9, A) is the 
Kullback-Leibler information number. Here, as the threshold 
is a time varying and decreasing function, the quantisation 
(0 is changed in the following way: if Off > 9* 



Yi 



k.l 



b\ if W k ,i 
b\ if W k ,i 
b\ if W kJ 



b\ if W k ,i G 



[g(kc),g(kc3A)), 
[g(kc3A),g(kc2A)), 
[g(kc2A),g(kcA)), 
[g(kcA),oo). 



(11) 



If 9 N < 9* we will transmit from {&", 0°, &°} under 
the same conditions. Here A is a tuning parameter and 
< 3 A < 1. The choice of d\ in (JHJ affects the performance 
of E[N\H Q ] and E[N\H X ] for the algorithm ©-(HQ. As 9 X 
increases, E[N\Hq] decreases and E[N\H{\ increases. 

The performance comparison of GLR-SPRT and GLR- 
CSPRT for the example in Section III C (with Z k - Af(0, 1)) 
is given in Table UTI Here A = 0.25. As the performance under 
Hi and H$ are different, we give the values under both. We 
can see that GLR-SPRT is always inferior to GLR-CSPRT. For 
Edd under H\, interestingly GLR-CSPRT have lesser values 
than that of SPRT-CSPRT for P FA > 0.02 (note that SPRT- 
CSPRT has complete knowledge of the SNRs), while under 
H it has higher value than SPRT-CSPRT. 



Hyp 




Pfa = 0.1 


P FA = 0.05 


Pfa = 0.01 


Hi 


SrKl-CSFKl 


1.615 


2.480 


4.28 


Hi 


GLR-SPRT 


1.597 


2.783 


5.286 


Hi 


GLR-CSPRT 


1.138 


2.221 


4.533 


H 


SPRT-CSPRT 


1.533 


2.334 


4.225 


H 


GLR-SPRT 


2.985 


4.257 


7.047 


Ho 


GLR-CSPRT 


2.424 


3.734 


5.72 



TABLE II 

Comparison among SPRT-CSPRT, GLR-SPRT and GLR-CSPRT for 

DIFFERENT SNR'S BETWEEN THE PRIMARY AND THE SECONDARY USERS 



The above scenario can also occur if the fading channel 
gain hi is not known to the Cognitive node /. Then under slow 
fading with Rayleigh distribution and using energy detector at 
the Cognitive Radios, / ,j ~ N(Q,a 2 ) and f ltl - N(9,(j 2 ) 
where 9 is random with exponential distribution and a 2 is 
the variance of noise. Here we use the GLR-CSPRT with the 
composite sequential hypothesis given in dH). The parameter 
9i is chosen as the median of the distribution of 9, such 
that P(9 > 9i) = 1/2. This seems a good choice for 9i 
to compromise between E[N\Ho] and E[N\Hi}. We use the 
example given in Section III C with Zk ~ A/"(0, 1) and 
9 ~ exp(l). Table [Ell] provides comparison of DualSPRT, 
GLRSPRT and GLRCSPRT. Notice that the comment given 
for the Edd in Table UTI is also valid here. 



Hyp 




Pfa = 0.1 


Pfa = 0.07 


Pfa = 0.04 


Hi 


DualSPRT 


1.74 


1.948 


2.728 


Hi 


GLR-SPRT 


1.62 


3.533 


9.624 


Hi 


GLR-CSPRT 


0.94 


1.004 


4.225 


Ho 


DualSPRT 


1.669 


1.891 


2.673 


Ho 


GLR-SPRT 


3.191 


3.849 


4.823 


Ho 


GLR-CSPRT 


2.615 


3.192 


4.237 



TABLE III 

Comparison among DualSPRT, GLR-SPRT and GLR-CSPRT with 

DLOW FADING BETWEEN THE PRIMARY AND THE SECONDARY USERS 

IV. Performance Analysis of SPRT-CSPRT 

Edd and Pfa analysis is same under Hi and Hq. Hence 
we provide analysis under Hi only. 

A. Pfa Analysis 

Let Po and Pi denote the probability measure under Hq 
and Hi respectively. Between each change of drift (which 
occurs due to the change in number of Cognitive Radios 
transmitting to the fusion node and due to the change in 
the value transmitted according to the quantisation rule 0) 
at the fusion center, under Hi, has a positive drift and 
behaves approximately like a normal random walk. also 
has a positive drift, but due to the min in its expression it 
will stay around zero and as the event of crossing negative 
threshold is rare becomes a reflected random walk between 
each drift change. The false alarm occurs when the reflected 
random walk crosses its threshold. Under Hi, let 



r l3 = inf{fc > 1 : F° > /?}. 



(12) 



We call Tp the first passage time at the fusion center. Let r 7) ; 
be the first passage time to threshold 7 by the I th node. Let tk 
be the k th order statistics of L i.i.d. random variables. Then 



Pfa at the fusion node, when H\ is the true hypothesis is 
given by, 

Pn 1 (False alarm) = Pn 1 (False alarm before t%) (13) 
+Pu 1 (False alarm between t\and tz) 
+ Pn 1 (False alarm between t 2 and t 3 ) + .... 

The main contribution to Pfa comes from the first 
term. 

P Hl (FA before h) 

CO 

k=l 

oo 

k=l 

In the following we compute Po\jp > x\rp < t\) and P[ti > 
k]. It is shown in £18] that, 

lim Pq{t/3 > x\t/3 < t±} = e~xjp(—\px),x > 0. (15) 

ft— >oo 

By finding solution to the integral equation obtained via 
renewal arguments [19], we can obtain the mean of first 
passage time, 773 (as done in (TJ, (2)). Let L(s) be the mean of 
T/3 with F ° = s and S k = log [ gi (Y k ) /g (Y k )} + D . From 
the renewal arguments, by conditioning on Sq = z: 

L(s) = F s (-s)L(0) + (16) 




L(s + z)dF s (z)dz + P[S> (3- s 



where F$ is the distribution of S k before the first transmission 
from the local nodes. By solving these equations numerically, 
we get X p = 1/L(0). 

Next we consider the distribution of t x . SPRT {W k ,i, k > 0} 
is a random walk at each secondary node I. We assume 
fo,i ~ 7V(0, of) and f u ~ Af(6 h af), where 6 t is the 
post change mean and af is the variance for I th Cognitive 
Radio. Let mean and variance of the drift of I th Cogni- 
tive Radio be St = E Hl [log (/1 (X u ) //„ (X k ,i))] , E? = 
Var Hl [log [/1 (X fe>i ) //o )]] respectively. We know Si > 
0. The time r 7i ; for W k j at each local node / to cross the 
threshold 7 satisfies i?[r 7j ;] ~ l/Si for large values of 7 
(needed for small Pfa)- Then by central limit theorem we 
can show that at each node / 

y z > 

Thus now (O equals 

00 L 

«2(l-e- A " fc )n(l-*r 7 ,(A)) 

k=l 1=1 

where $ r , is the Cumulative Distribution Function of r 7i j, 
obtained from the Gaussian approximation 1171 . 

Table [IV] provides comparison of Pfa via simulation and 
analysis. 



7 


/3 




-fV.4 Anal. 


E DD Sim. 




15 


30 


0.0072 


0.0065 


33.1585 


31.7624 


12 


27 


0.00675 


0.00613 


26.8036 


24.9853 


14 


26 


0.01675 


0.01624 


30.0817 


29.1322 



TABLE IV 

SPRT-CSPRT FOR DIFFERENT SNR'S BETWEEN THE PRIMARY AND THE 
SECONDARY USERS: COMPARISON OF E DD AND P FA OBTAINED VIA 
ANALYSIS AND SIMULATION. UNDER Hi 



B. Edd Analysis 

In this section we compute Eod theoretically, ti, i th order 
statistics of L random variables r 7 ./. 1 < I < L, is the first 
time at which i local nodes are transmitting. Mean of ti can 
be find out from the method explained in [3], for the finding 
k th central moment of non i.i.d. L th order statistics. 

Between ti and U + i the drift at the fusion center is not 
necessarily constant because there are four thresholds (each 
corresponds to different quantizations) at the secondary node. 
The transmitted value changes after crossing each threshold, 
61 — > i>2 • • • — > 64. Let 1 < j < 3 be the time points at 
which a node changes the transmitting values from bj to bj + i 
between ti and £j+i. We assume that the probability of false 
alarm at the local nodes, Pf a is very small. Also with a high 
probability the secondary node with the lowest mean in dT7| i 
will transmit first, the node with the second lowest mean will 
transmit second and so on. In the following we will make 
computations under this approximations. The time difference 
between ip* and tpP +x transmission can be calculated if we 
take the second assumption (=Ai/(5;)- We know E[U] for every 
i from an argument given earlier. Suppose I th node transmits 
at tf instant and if E[U] + A 1 /S l < E[t l+1 ] then E[t\] = 
E[U] + A-l/Si. Similarly if E[t\] + A 1 /S l < E[t l+1 ] then 
E[tf\ — E[tj]+Ai/Si and so on. Let us represent the sequence 
t = {ti,t\,tf,ti,t2, •■•,£5} (entry only for existing ones by 
the above criteria) by T = {T U T 2 ,T 3 , ...}. 

Let fi k be the mean drift at the fusion center between T k and 
T k+ i. Thus Tfc's are the transition epochs at which the fusion 
center drift changes from u k -i to fj, k . Also let F k = E[Fx k -i] 
be the mean value of F k just before the transition epoch T k . 
With the assumption of the very low Pf a at the local nodes and 
from the knowledge of the sequence t we can easily calculate 
li k for each T k . Similarly F k+1 = F k +fi k (E[T k+1 } - E[T k ]). 
Then, 

E DD w E[T 3 ] + ^— ^ (18) 

where 

/3 - Fi 

j = min{i : m > and < S[Tj+i] — S[Tj]}. 

Mi 

Table [TV] provides the simulation and corresponding analysis 
values. We used the same set-up as in Section III C (with 

Z k ~ 7V(0, 1)) 

V. Conclusion 

We consider the problem of cooperative spectrum sensing 
in this paper. We provide improved algorithms SPRT-CSPRT 
and DualCSPRT over a recent algorithm DualSPRT. We show 
that these algorithms can provide significant improvements. 



We provide theoretical analysis of SPRT-CSPRT and compare [28] Q. Zou, S. Zheng, A. Sayed, "Cooperative Sensing via Sequential 
to Simulations. We further extend these algorithms to cover Detection", to appear in IEEE Trans. Signal Pwces., 2010 

the case of unknown SNR and channel fading. 
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